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We discuss how to analyze the optical conductivity and Raman spectra of multi-orbital systems 
using the velocity and the Raman vertices in a similar way Raman vertices were used to disentangle 
nodal and antinodal regions in cuprates. We apply this method to iron superconductors in the 
magnetic and non-magnetic states, studied at the mean field level. We find that the anisotropy in 
the optical conductivity at low frequencies reflects the difference between the magnetic gaps at the 
X and Y electron pockets. Both gaps are sampled by Raman spectroscopy. We also show that the 
Drude weight anisotropy in the magnetic state is sensitive to small changes in the lattice structure. 

PACS numbers: 75.10.Jm, 75.10.Lp, 75.30.Ds 



I. INTRODUCTION 

Raman and optical conductivity are very useful tech- 
niques to analyze the electronic properties of strongly 
correlated electron systemsPE] Optical conductivity ex- 
periments have provided very valuable information on the 
reorganization of the spectral weight and the opening of 
gaps in many materials.^ In cuprates, the use of different 
polarizations in Raman scattering has allowed the disen- 
tanglement of the differen t ph ysics of the nodal and the 
antinodal electronic states!^ 

The multiband character of iron superconductors com- 
plicates the analysis of their RamarP^ anc } optical 
conductivity^^!] spectra. The five iron 3d orbitals are 
required for a minimal model to describe these materi- 
als. Different interband transitions involve similar ener- 
gies and contribute to the spectra in the same fr equ ency 
range. Moreover they start at very small energies^ mak- 
ing it difficult to separate their contribution from the 
Drude peak in optical conductivity experiments! 13 l 2 ° l 

The difficulties in the interpretation of the spectra are 
more pronounced in the magnetic state. When entering 
the magnetic state the optical conductivit y is suppressed 
at low frequencies and new peaks appear P^ 2 ! * 14 ! 16 ! 18 * 21 ! 
The conductivities along the antiferromagnetic a xx {uj) 
and ferromagnetic a yv (w) direct ions show different inten- 
sity and peak frequencies! 16 * 21 ^ The modification of the 
spectrum and the anisotropy in the magnetic state are 
visible up to 17000 cm -1 , Ref. 21 . The in-plane re- 
sistivity in the magnetic state is also an isotropic . The 
origin of these anisotropies is not clear yet! 16 * 21 *^^ The 
Raman spectrum in the magnetic state shows signatures 
and peaks at energies similar to those found in optical 
conductivity and it has been interpreted in terms of two 
kinds of electronic transitions: a high energy transition 
between folded anti-crossed spin density wave bands and 
a lower energy transition which involves a folded and a 
non-folded bandP see Fig. [l] It would be desirable to ad- 
dress the orbital degree of freedom in the interpretation 
of the optical conductivity and Raman spectroscopies. 

In this paper we discuss how to analyze the optical 
conductivity and Raman spectrum of multi-orbital sys- 
tems using the velocity and the Raman vertices. These 



vertices depend on the symmetry of the orbitals involved 
in the interband transitions but do not simply follow the 
symmetry rules for atomic transitions. They change in 
k-space, see Fig. [2j reflecting the underlying lattice and 
can be used to obtain information of the orbitals and 
the regions of the Brillouin zone (BZ) which contribute 
to the spectrum in a similar way Raman vertices were 
used to disentangle the nodal and antinodal regions in 
cuprates. We apply this method to iron superconductors 
in the magnetic and non-magnetic states. The informa- 
tion obtained from optical conductivity and Raman spec- 
troscopy are complementary and allow the exploration of 
all the BZ. 

We find that for magnetic moments comparable to the 
experimental ones the different frequencies at which a xx 
and (J vy peak reflect the magnetic gaps at Y and X elec- 
tron pockets respectively (in the one-Fe unit cell). For 
some reconstructed band structures any of these gaps can 
open below the Fermi level. In this case the interband 
transition is not allowed affecting the shape of the spec- 
trum. We also show that the Drude weight anisotropy is 
sensitive to small changes in the lattice structure. 

The article is organized as follows: in Section [n] we 
give the expressions for the optical conductivity and the 
Raman scattering in multi-orbital systems and introduce 
the respective vertices. Section [TTT] focuses on the Raman 
and velocity vertices in the case of iron superconductors. 
In Section |IV| we use the vertices to discuss the optical 
conductivity and Raman spectra of iron superconductors 
in the mean- field magnetic and non-magnetic states. The 
fingerprints in the spectra of the crossover to the orbital 
differentiation regime which appear in our mean field cal- 
culations are analyzed. Section [V] is dedicated to the 
Drude weight in-plane anisotropy. We end with a discus- 
sion of our results and a comparison to experiments in 
Section EH 
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FIG. 1. (Color online) Illustrations of the types of transitions 
encountered: (a) standard interband transition; (b) transi- 
tions involving a non-folded and a folded band; (c) transitions 
involving two folded anticrossing bands. 



man response. V is the volume, the Fourier com- 
ponent of the electromagnetic field, and e a the photon 
polarization with energy w q while e Q is its complex con- 
jugate. 



A. Optical conductivity 

Assuming a single component of the vector potential 
A° and expanding to linear order, the longitudinal cur- 
rent is given by 



II. MODEL AND METHOD 

We consider a multiorbital system with a tight binding 
Hamiltonian 



Hn = 



e Ml/ (k)c kAicr Ck„ 



(1) 

i and j label the lattice sites connected by the hopping 
terms, /i and v the orbitals, and a the spin. e^(k) is the 
tight binding in k-space. 

If interactions are local, as considered through all 
the paper, the coupling between the electrons and the 
electromagnetic field can be introduced via the Peierls 
substitutiorP2lIini in the hopping if/ -> ^.V JV^M ^ 

I = l x x + l y y links the neighbors i and j connected by 
hopping. Here and in the rest of the paper natural units 
are used and we take e = c = H = 1. Assuming that 
the vector potential A(r) varies more slowly than the 

lattice length scale one can approximate f? +l A(r) ■ df ~ 
A(i+i/2) • I. For small fields we can expand the exponen- 
tial to second order e lAl ~ 1+iA-l— ^(A-l) 2 . Under these 
approximations and introducing Fourier components, the 
Hamiltonian in the presence of an electromagnetic field 
is 



H(A)=H(A = a)-J2Y,% A -i 

a q 

- - V V T aP A a At 



q 2 ' 



with a, f3 = x,y and 
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dk a dkp c k+§^ Ck -§ i 
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M + ''q+q'^-q' 

q' 



(5) 



The contribution of the paramagnetic current to the 
expected value of the current is given by the Kubo for- 
mula. At zero temperature 
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where \<po) is the ground state. We are interested in 
the response to a homogeneous q = electric field 
A^ =0 (cj) = Eq =0 (cj) / (iw — S) with 5 a small parameter, 
and in particular in a' aa (uj), the real part of the opti- 



•^q=o = o-aa(w)^" (w). After some algebra,™Sl 



cal conductivity. (t,„,u-) = o4 a (w) + i(7^ a (w^d< > liiu > <l as 
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with the Drude weight given by 
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Eq and 2£ m are the energies of the ground state |</>o) and 
the excited states \<f> m ) respectively. The first and second 
terms in Eq. Q originate respectively in the diamagnetic 
and paramagnetic contributions. Eq. ([7| fulfills the op- 
tical sum rule^l 



a' aa {uj)duj = 7t(0 o 



rpOLOL 

~ J q=o| 



(9) 



that becomes the kinetic energy when hopping is re- 
stricted to first nearest neighbors. 

Within the mean-field level used below^H the Hamilto- 
nian becomes biquadratic in fermionic operators. There- 
fore, the eigcnstates \4> m ) can be given in terms of single 
particle bands and Eqs. ^ and Q can be written 



Aa 



(4) 



jq is the paramagnetic current and gives rise to the 
diamagnetic current and enters in the non-resonant Ra- 



o> aa (w)=D a 6(u>) + ^ 



x0(e„,(k))0(-e„(k))«S(u; - e„,(k) + e„(k))(10) 
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FIG. 2. 



(Color online) Leading k-dependence of the velocity (optical conductivity) and Raman squared vertices \v* 
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and 



defined in Eqs. (24h(26l for iron superconductors for each pair of orbitals involved in an interband 



transition in the unfolded 1-Fe Brillouin zone corresponding to the tight-binding model in Ref. [38]. x and y directions point 
along the Fe-bonds. B\ g and B2 g are defined in this Brillouin zone and axes. The same assignment is used in the non-magnetic 
and the magnetic state. In the case of intraorbital transitions B\ g probe the electron pockets. Note that a different convention 
has been used in some Raman articles which define the polarizations according to the tetragonal FeAs Brillouin zone. Darker 
regions in each figure correspond to larger values of each squared vertex. Combination of optical conductivity and Raman 
scattering spectroscopy allow to probe selectively inter-orbital transitions through the whole Brillouin zone. Intraorbital B\ g 
squared vertex depends on k as {Acosk x — Bcosk y ) 2 with the ratio A/B dependent on the orbital. A = B valid for xy, 
x 2 — y 2 and 3z 2 — r 2 has been considered in the figure. The k-dependence changes slightly for zx, yz, but it is maximal in 
approximately the same regions of the Brillouin zone . The pairs (±,±) refer to the product of the parities with respect to x 
and y reflections for each pair of orbitals, see Sec. Ill for an illustrative example. The vertices for the (+, +) 3z 2 ~r 2 -t^r x 2 — y 2 
transition in B-z g symmetry and for the (— , — ) transitions in B\ g symmetry vanish. Note that the hoppings between xy and 
x 2 — y 2 to nearest and next nearest neighbors are zero, so the corresponding vertices cancel within the model considered here. 
A finite small contribution could arise if hoppings to farther away neighbors were considered. 
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A,= -7r^C(k)0(-e„(k)) 



2^ v liZJMf 
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kn^n' x ' v 



0(e n ,(k))0(-e n (k)),(ll) 



where e„(k) are the band energies, 0(e n (k)) the Heaviside 
step function and 

O k ) = ]T ^^< n (k)a„ n (k) , (12) 



#„(k)= - E %^<n'( k Kn( k ) , (13) 

with a M „ (k) the rotation matrix between the orbital and 
the band basis cj^ = E„ ajn(k)4ncr- 

B. Raman response 

The electronic Raman scattering measures the total 
cross section of the inelastic scattering of electrons. It is 
proportional to the transition rate of scattering an inci- 
dent (qj, LOi, ej) photon into an outgoing (q s , lu s , e s ) state, 
where qi jS , uii yS and ei yS label the momentum, frequency 
and polarization of the incident and scattered photons. 
The transi tion rate can be obtained following the Fermi 
golden rultP2l 



^e-^ E '\{F\M\I)\ 2 5{E F -Ej 



UJ 



(14) 



I.F 



with Z the partition function and /3b = 1/KbT, being 
Kb the Boltzmann factor and T the temperature. E] 
and Ep are the initial and final energies of the many 
electron system, u> = u s — W{ is the transferred energy 
and M is the effective light scattering between the initial 
and final state. Neglecting resonant processes^ 



(F\M\I)=J2^ S (F\ 



q I 



i) 



(15) 



a/3 



with q = q s q^. For the energies involved in the Raman 
scattering q << k^?. In the following we take q = 0. 

Instead of using T^£ and arbitrary a and j3 polariza- 
tions, it is convenient to decompose this matrix clement 
into basis functions of the irreducible point group of the 
lattice^ according to the polarization of the incident and 
scattered light A= Bi g , B2 g , etc and use an effective po- 
larization dependent density matrix 



P X =Y, 7^(k)c 



k.cr [if 



k/icr 



In particular, for Bi g and B 2g polarizations 

^/u^^VW d%„(k) 
7m" 1 k: 



dkl 

7 £^(k)= d 2 e^(k) 



(16) 

(17) 
(18) 



At zero temperature the Raman spectrum for A polariza- 
tion becomes: 

S X (fj) (xJ2\( ( t ) rn\p X |0 O )| 2 6{E m - E - u) . 
m 

If, as in the previous subsection, the eigenstates can be 
expressed in terms of single-particle bands with energies 
e»(k) 

Sx(w) = J2 hn'n(V\ 2 e(- e „(k))6(e n ,(k)) 



xS(e n /(k) - e„(k) - uj) , 

Jn'n ( k )= J2 Tp* ( k Kn> ( k K™ ( k ) ■ 



(19) 



S\(u>) is related^) to x A (w), the imaginary part of the 
effective polarization density correlation function Xx(w), 
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discussed below by: 

Xx{uj) = 7T [Sx(w) - S\(-u)] 



<[/(*), /(0)]>, (20) 



(21) 



At zero temperature and zero scattering rate S\(oj) oc 
X A (w). In our calculations we broaden the delta functions 
with a small scattering rate T = 20 meV. With such T, 
the proportionality relation between Xa( w ) an< ^ ^a( w = 

0) only fails at u ~ 0. Xa( w = °) = wllilc S ^ = °) 
acquires a small finite value. However, the qualitative 
features of the spectra are not affected. 



C. Raman and velocity vertices 



According to Eqs. (10 1, (131 and (19) the optical con- 
ductivity and Raman spectrum are given in terms of the 

2 I \ I 2 

quantities |i" n /(k)| and 7„ n /(k) . Once the square 
is expanded in orbitals the optical conductivity involve 
terms 



^(k)Se* v ,(k) . 



^ dk a 



dk a 



"fin 



(k)ay n (k)«w (k)a*,„, (k) . 



Similar terms appear in the calculation of the Raman 
spectrum. The transition probability between n and n' 
has a very complex fc-structure which depends on the 
interactions through the a^ n operators, making very dif- 
ficult the analysis of the spectra. 

In order to get a simple picture of the transitions, one 
can neglect the crossed terms with fj, ^ // or v ^ v' and 
focus on the terms with /i = // and v = v 1 . Within this 
approximation the probabilities for an interband tran- 
sition between bands n and n' in optical con duct ivity 



\jn n '\ an d Raman spectroscopy 7„„/ in Eqs. (10) and 
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FIG. 3. (Color online) Upper figures: Energy bands in the non-magnetic (a) and (n, 0) antiferromagnetic state (c) for U = 
Jh = and U = 1.6 eV and Jh = 0.25(7 respectively. Linewidths and colors reflect the orbital content yz=red, za;=green, 
xy=bhxe, 3z 2 — r 2 =orange and x 2 — y 2 =black. Note that not all the orbital content is visible due to overlapping curves. Lower 
figures: Zoom of the band structures in the upper figures in the energy region close to the Fermi level showing the optical and 
Raman transitions at low energy discussed in the text. 



^ are 

|^„,(k)| 2 ~ J2 IC( k )f K«(k)| 2 |<w(k)| 2 , (22) 



Inn 



'( k )f ~ E W,(k)| 2 1^"^! 2 l°-'( k )| 2 ■ ( 23 ) 



component fi in band n to the orbital component v 
in band n' in a Raman or optical conductivity experi- 
ment are respectively weighted by the squared vertices 
|7^„(k)| and ^^(k)] . These vertices depend on k 
through the underlying lattice via e M „ (k) . While the ver- 



tices \v" (k)| and 



7M" 



(k) 



depend on the symmetry 



Here |a Mn (k)| and |a l/ „/(k)| give the spectral weight of 
orbitals fj, and v in the bands n and n' involved in the 
transition. For each transition between bands n and n' , 
this approximation only keeps the weights corresponding 

to fj, orbital-density on band n |a M „| 2 to v orbital-density e^(k) in the orbital basis helps to identify which transi 
on band n |aw| 2 , summing to all pairs of orbitals. The 
velocity | ^^(k) | and Raman Vy^v (k) | squared vertices 
are given by 



flu V 

of the orbitals, they cannot be deduced from simple ar- 
guments involving atomic optical transitions, valid only 

at r. 

Therefore, knowledge of the tight binding dispersion 



k„(k)| 



7^ l9 (k) 



7m. 29 (k) 



de M „(k) 



dk a 



c»V( k ) d 2 e M „(k) 



dk 2 x 
d%„(k) 



dkl 



(24) 



(25) 



(26) 



)1U\ 

tions contribute to the optical conductivity and Raman 
spectrum. The band structure of the interacting Hamil- 
tonian in the normal or ordered state is also required for 
the interpretation of the spectrum via the orbital spectral 
weight |a„„| 2 and the band energies e n (k). 
We emphasize that the approximation in Eqs 



22p6~ 



where a — x,y and we have focused on the A = B\ g , B2 g 
Raman polarizations. 

Besides the corresponding spectral weights la^j 2 and 
|cw| 2 , the probabilities of a transition from the orbital 



is used only in the interpretation of the spectra that is 
calculated using the full expressions in Eqs. [T0| and [T9| 



III. VERTICES IN IRON SUPERCONDUCTORS 

Previous expressions are valid for a multiorbital sys- 
tem with local interactions. In the rest of the paper we 
focus on the case of iron superconductors. Unlike oth- 
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erwise indicated we consider the five-orbital tight bind- 
ing model introduced by the authors in Ref. [38], where 
the five orbitals refer to the 3d iron orbitals yz, zx, xy, 
3z 2 —r 2 , x 2 — y 2 . The model and the orbital directions are 
defined in the one-Fe unit cell, with x and y along the 
Fe-Fe bonds. Hopping among the orbitals is restricted 
to first and second neighbors and includes both indirect 
hopping mediated by As as well as direct hopping be- 
tween Fe orbitals. Indirect hopping depends on the angle 
&Fe-As formed by the Fe-As bond and the Fe plane. The 
tight-binding model used allows the modification of this 
angle. Except otherwise indicated the results are given 
for a regular tetrahedron with ctFe-As = 35.3° with a 
non-magnetic band structure as plotted in Fig. [3] (a). 
Similar results are expected for other five-orbital models 
discussed in the literature, whenever defined using the 
same orbitals and unit cell. 

The five orbitals result in 15 different squared vertices 
(5 intraorbital and 10 interorbital) for each of the spec- 
troscopies discussed in this article: a' xx , a' yy , XB lg anc ^ 
■ Keeping only the leading k-dependence, one can 
group the squared vertices according to the symmetry 
of the product of the orbital wavefunctions involved in 
the transition with respect to x and y reflections. For 
example, yz orbital is even (odd) with respect to x (y) 
while xy is odd with respect to both x and y. There- 
fore, the product yz ■ xy is odd (even) with respect 
to x (y), summarized in (—,+). The pairs yz,xy and 
zx, 3z 2 — r 2 have the same product parity (—,+). For 
both terms the derivative involved in the vertex of a y 
is de^u/dky — 2icosk y (tJ ll , — 2£ MjJ/ cos fc^) with t and 
t the hoppings for the corresponding orbitals to nearest 
and next nearest neighbors respectively!^ The leading 
/c-dependence of the squared vertex for a y is cos 2 k y for 
the orbital pairs with product parity (—,+). 

The leading k-dependence of the vertices is shown in 
Fig. [2] where the darkest color emphasizes the region 
of the BZ with largest square vertex and the product 
symmetry is given in parenthesis (±, ±) with + and — 
referring respectively to even and odd. 

In Fig. [2] the vertices for the (— ,— ) transitions in 
Big symmetry vanish. The (+, +) symmetry group in- 
cludes both the transition 3z 2 — r 2 «-» x 2 — y 2 and 
the intraorbital ones. The optical conductivity veloc- 
ity vertices of these transitions have the same leading 
k-dependence but there are some differences in the Ra- 
man case. While the squared B 2g vertex vanishes for the 
3z 2 — r 2 f-> x 2 — y 2 transition it depends as (sin k x sin k v ) 2 
for intraorbital transitions. The B\ g squared vertex goes 
like (Acosk x — Bcosk y ) 2 with A and B orbital depen- 
dent. A = —B for the 3z 2 — r 2 ■<-» x 2 — y 2 transition while 
A = B for intraorbital transitions involving xy, x 2 — y 2 
or 3z 2 — r 2 . This results in a different k-dependence of 
the squared vertices with maximum and minimum values 
in different regions in k-space. For zx and yz orbitals 
A 7^ B but the k-dependence is similar to that shown for 
A = B. 

The x 2 — y 2 intraorbital Raman squared vertices de- 
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FIG. 4. (Color online) Main figure: (n, 0) magnetic phase of 
the undoped (6 electrons in 5 d orbitals) system as a function 
of the interactions U and Jh/U. The white area is the non- 
magnetic region. The blue and red areas are magnetic with 
a high moment (parallel orbital moments) and a low moment 
(antiparallel orbital moments) state respectively, studied in 
Ref. [41] . The dashed line separates the itinerant and the 
strongly orbital differentiated magnetic regions, see text and 
Ref. |44j for discussion. Inset: Doping dependence of U mag 
(U*), the interactions at which antiferromagnetism (strong 
orbital differentiation) appear for Jh = 0.25(7, reproduced 
from Ref. |H]. 

pend on k as (cos k x — cos k y ) 2 for B\ g and (sm.k x sin k y ) 2 
for B 2g - This dependence was widely used in the analysis 
of the cuprates Raman spectrum and allowed to separate 
the physics of the antinodal region sampled by B\ g from 
that of the nodal region measured by B 2g . In a similar 
way as done in the cuprates for the x 2 — y 2 intraorbital 
transitions, the contributions to the spectrum of a tran- 
sition between two bands weighted by two given orbitals 
will be larger when the corresponding vertex is large in 
the region of k-space which satisfies the energy conser- 
vation and this can be used to analyze the experimental 
and calculated spectrum. The region of the BZ weighted 
for transitions between two orbitals is different for the 
optical conductivity and Raman vertices, being possible 
to cover practically all the BZ for every transition by 
studying a' xx , u' yy , x' Blg and Xb 2s - Note that due to the 
multiorbital character B\ g [B 2g ] can sample different re- 
gions of the BZ besides the antinodal (nodal) region. 



IV. OPTICAL CONDUCTIVITY AND RAMAN 
SPECTRUM OF IRON SUPERCONDUCTORS 

A. Phase diagram and band reconstruction 

To discuss the spectrum of iron superconductors in the 
magnetic and non-magnetic states, we consider a five 
orbital interacting Hamiltonian with local interactions. 
Assuming rotational invariance, the interactions can be 
written in terms of two parameters: the intraorbital in- 
teraction U and the Hund's coupling Jh, see Ref. [H] for 
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FIG. 5. (Color online) Band structure in the Y-Y-X directions 
in the antiferromagnetic state corresponding to Jh = 0.25(7 
with (a) U = 1.6 eV and (b) U — 1.8 eV, in the itinerant 
and orbital differentiated regions respectively. Linewidths and 
colors reflect the orbital content yz=red, za;=green, xy=b\ue, 
3z 2 — r 2 =orange and x 2 — y 2 =black. The circles highlight 
the shift of the yz orbital to higher energies, driving yz to a 
half-filled gapped state, and x 2 — y 2 to lower energies. The 
arrows mark the involved transitions. See text for discussion. 



details. Except otherwise indicated we assume electron 
filling n = 6, as in undoped materials, and use our tight 
binding model in Ref. [38] with squared vertices as given 
in the previous section. The Hamiltonian is treated at 
the mean field level. 

The mean field (it, 0) magnetic phase diagram as a 
function of the i nterac tions U and Jh/U has been dis- 
cussed previousl y 41 * 44 ! and is reproduced in Fig. |] for 
clarity, see also Refs. [IS] and The white area is 

the non-magnetic state and the red area corresponds to 
a low moment state with antiparallel orbital moments 
violating Hund's rulePD The blue area corresponds to a 
magnetic state with parallel orbital moments. In this 
phase there is a crossover from an itinerant to a strong 
orbital differentiation regime represented in Fig. [4] by a 
white dashed line. In the strong orbital differentiation 
regime xy and yz are half-filled gapped states while zx, 
3z 2 — r 2 and x 2 — y 2 are itinerant with a finite density of 
states at the Fermi levelP^ 

With hole doping the system becomes more 
correlatecPESEo] as the average orbital filling ap- 
proaches half-filling. As shown in the inset of Fig. [4] 
the interaction U* at which the system enters into the 
orbital differentiated regime decreases with hole doping 
while U magi the interaction at which antiferromagnetism 
appears, is non-monotonous™ 

Fig. [3] (c) shows the reconstructed bands in the (n, 0) 
magnetic itinerant regime for U = 1.6 eV and Jh = 
0.25U corresponding to a magnetic moment m = 0.91u,g. 
The band reconstruction involves band foldings not only 
at the Fermi surface but also far from it. Gaps with 



different values open at different points in the Brillouin 
zone in the bands below, at and above the Fermi level. 
Breaking of the zx — yz degeneracy is observed. Along 
r — X new bands cross the Fermi level forming V-shaped 
pockets referred to in the literature as Dirac pockets. 

When moving towards the orbital differentiated phase 
with increasing interactions or decreasing doping, the 
magnetic moment increases and the band reconstruction 
becomes more complex. Bands become flatter and shift 
away from the Fermi level, being this effect more accused 
for the yz and xy orbitals. Spectral weight of these two 
orbitals shifts partially to higher energies above the Fermi 
level as these orbitals become half filled. As shown in 
Fig. [5] (a) for U = 1.6 eV and J H = 0.25Z7, in the itin- 
erant regime the yz band at Y slightly above the Fermi 
level (red) lies below a band with dominant x 2 — y 2 or- 
bital content (black). On the contrary for U = 1.8 eV, 
Fig. [5] (b) , in the orbital differentiation regime the order 
of these two bands is reversed, driving yz to half-filling. 
This feature seems to be a fingerprint of the orbital differ- 
entiated regime, at least in all the cases analyzed for our 
tight binding modeP^ and the one by Graser et al [51] . 



B. Optical conductivity 

The optical conductivity is shown in Fig.[6]for the non- 
magnetic state corresponding to U = Jh = (dotted 
lines) and for the magnetic state for U = 1.6 eV and 
J H = 0.2517 (solid lines) using our modelP Fig. ga), 
and the one by Graser et alp' Fig. J^b). A scattering 
rate of 20 meV is assumed in the calculations. 

Due to the tetragonal symmetry in the non-magnetic 
state (t' xx (w) = a' (w). Beyond the Drude peak, at low 
frequencies, interband transitions are clearly visible start- 
ing from 0.35 eV and contribute to the conductivity up 
to energies above 5 eV. In the following we focus on the 
low energy transitions. As marked in Fig. [3jb) and fol- 
lowing the squared vertices in Fig.[2j the peaks in o~ xx (u>) 
and o~ vy (uj) around 0.4 eV originate in transitions at Y 
involving 3z 2 — r 2 — > yz and 3z 2 — r 2 — > zx respectively. 
Similarly, the shoulders at 0.45 and 0.7 eV in a xx {oj) come 
from interband transitions involving zx — > xy along Y — Y 
and yz — > x 2 — y 2 close to Y, but not at Y. Symmetry re- 
lated transitions yz — > xy along X — Y and zx — > x 2 — y 2 
close to r, but not at Y, give an identical contribution 
to 

<J vy{ L0 )- While not easily identifiable in the spectrum, 
interband transitions between the two hole pockets at Y 
contribute at frequencies within the Drude peak. This 
interband transition is very narrow, as can be seen in the 
inset of Fig. [6]ja) where the Drude peak has been sub- 
tracted from the optical conductivity. This transition is 
allowed by the finite e g content of the yz j zx dominated 
hole bands P^l 

In the (tt,0) magnetic state a' xx (to) ^ a yy (uj), as ex- 
pected from the tetragonal symmetry breaking. The two 
curves cross many times through all the frequency range, 
namely, the sign of the optical conductivity anisotropy is 
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FIG. 6. (Color online) Optical conductivity for two differ- 
ent tight-binding models. A logarithmic scale for to is used 
to highlight the low energy region. A scattering rate of 20 
meV is used, (a) Optical conductivity for our tight binding 
modelP^ The dotted line corresponds to the non-magnetic 
state for U — Jh = eV. Due to the tetragonal symmetry 
a' xx {u)) = a' yy (u}). In solid lines, the optical conductivity in 
the antiferromagnetic x (black) and ferromagnetic y (red) di- 
rections corresponding to U = 1.6 eV and Jh = 0.25(7 in 
the (tt, 0) antiferrromagnetic state. The tetragonal symme- 
try is broken and a' xx (uj) ^ a yy (uj). There is a strong sup- 
pression of the Drude weight with magnetism: a' xxyy (uj = 
0,U = 0) = 12.0876, o' xx (u = 0, U = 1.6eV) = 2.2855 and 



0,(7 



1.6 eV) = 1.9874. Inset: Optical conduc- 



tivity at low energy in the non-magnetic state excluding the 
Drude peak. A peak corresponding to an interband transition 
at very low energies (lj ~ 50 meV) shows up. (b) Same as 
(a) for the tight binding model of Graser et al Ref. [51] with 
U — Jh = in the non-magnetic state and U = 1.42 eV 
and Jh = 0.25(7 in the magnetic state. See text for discus- 
sion. The Drude weights for the varios cases are: (j' xx yy (uj = 
Q,U = 0) = 6.679, a' xx {uj = 0, U = 1.42 eV) = 2.778 and 
a' yy (uj = 0,(7 = 1.42 eV) = 1.2062. 



frequency dependent. In the itinerant regime, for mag- 
netic moments ~ 0.9/is similar to the experimental ones, 
the main effect of magnetism is the partial suppression 
of the Drude peak, which becomes anisotropic, and the 
appearance of a magnetic peak, see Fig. [6] (a). This peak 
originates in transitions across a magnetic gap between 
two anticrossed bands. The position of this peak depends 
on the magnetic moment and is different for a xx (uj) and 
Oyy (bj) because, as it can be seen in Fig. [2j each conduc- 
tivity samples a different region in k-space. 

The magnetic peak in a yy {u) originates in the xy — > yz 
transition at the electron hole pocket at X along the 
X — M direction and the folded yz —> xy transition at 
the yz hole pocket at T along T — Y, see Fig. [3] (d). On 
the other hand, the peak in a xx (uj) is due to a zx — > xy 
transition close to Y along Y — M . A transition with the 
same energy happens at the hole pocket in M. Thus the 
peak in a yy (ui) measures the gap at the electron pocket 
at X and the hole pocket at T while the one in ct xx (uj) 
measures the gap at the electron pocket at Y which is 
the same as at the hole pocket at M. 

For the values reported in Fig. [6]ja) the magnetic peak 
in ct xx (lu) appears at larger frequencies than that in 
a yy {uj). The gap at the electron pocket at Y depends 
on the size and height of the hole pocket at M which 
is ve ry sen sitive to the position of the As in the FeAs 
layer I2SI521 pjg Jsjb) shows the optical conductivity corre- 
sponding to Graser et al tight-binding model in Ref. [5T] 
for U = 1.42 eV and J H = 0.25U. The magnetic mo- 
ment is m = 0.91/is, as for the values used with our 
tight-binding model in Fig. |6ja). However the positions 
of the magnetic peaks are reversed, with the magnetic 
peak on a xx at a lower frequency than that on a yy in the 
Graser et al case. The electronic bands corresponding to 
this model are shown in Fig. [7J The xy hole pocket at 
M touches the Fermi level in the non-magnetic bands. 
Even though the bands are similar in the paramagnetic 
state, compare Fig. [3][a) and Fig. [7^ a) , the reconstructed 
bands in the Y — M direction close to the Fermi level in 
Fig.JTTd) are very different from the ones in our model 
Fig.|3jd). The gaps at the electron pockets differ in both 
models even though we have chosen parameters such that 
the magnetic moment is the same. 

In the two cases discussed here the magnetic transi- 
tions at both electron pockets are allowed. However, if 
the minimum of any of the upper folded bands lies be- 
low the Fermi level the corresponding transition would be 
forbidden and the spectrum strongly modified. This hap- 
pens in Ref. |35j where a small electron pocket is formed 
close to X in the X — M directions (and correspondingly 
close to T in the V ~Y direction) resulting in the sup- 
pression of the a yy magnetic peak observed in our model 
and highlighted in Fig. [3^d). It can also happen in other 
models close to Y if the hole pocket at M is below the 
Fermi level in the non-interacting bands. 

The magnetic peak is the main but not the only sig- 
nature of magnetism in the optical conductivity at low 
frequencies. There are smaller peaks and plateau like 
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structures which originate in transitions close to T (and 
equivalently clo se to X) between a folded band and a 
non- folded bandf^H see Fig. [ljb). Whether these tran- 
sitions are allowed depends on the starting tight-binding 
model and the magnetic moment. For example, in our 
model, see Fig. ^d) , the gap in the zx hole pocket along 
r — X opens below the Fermi level at a value of mo- 
mentum k x for which the upper yz band is occupied and 
the transitions involving these bands at this k-point are 
forbidden. On the contrary, in Graser et al model, see 
Fig. [7^d) , these transitions are allowed and produce the 
peaks at 0.21 and 0.55 eV in <j yy (u>). The small contri- 
bution to a xx (uj) below the magnetic peak comes, on the 
other hand, from zx — > xy transitions along T — Y. 

With increasing magnetic moment, in the orbital dif- 
ferentiation region, the shape of the spectrum changes, 
see Fig. |8j and the spectral weight shifts to higher 
energies P^The interband transition between the 3z 2 — r 2 
band below the Fermi level at T and the yz band, strongly 
affected by the orbital differentiation in Fig. [5] is active 
for <7 xx (uj). It produces a step like feature in the spec- 
trum. However because the spectrum is very sensitive to 
parameters it is not possible to signal an easily identifi- 
able fingerprint of this transition in the spectrum. 



C. Raman spectra 

The B\g and Bi g Raman spectra in the non-magnetic 
state are plotted in Fig. [9^ a) . Both symmetries display a 
peak at very small energies, however the nature of these 
two peaks is different. In £? lg it corresponds to intra- 
orbital transitions at the electron pockets related to the 
Drude contribution to the optical conductivity, and it is 
allowed by the finite scattering rate. In Bi g it originates 
in the interband transition between the zx and yz hole 
pockets close to T. The rest of the spectrum comes from 
interband transitions. The strong peak in Bi g around 
0.7 eV corresponds to a transition between the 3z 2 — r 2 
and the xy hole pockets at M . The energy of these hole 
bands is very sensitive to the As height in the FeAs layer 
and the transition could be absent if both bands lie below 
the Fermi level. The peak at 0.4 eV in B ig originates in 
an interband transition at T between the 3z 2 — r 2 hole 
band and the hole pockets with non-negligible x 2 — y 2 
content. For higher energies there is a large bump in 
Big. It starts with a step like feature coming from the 
transition at V between the 3z 2 — r 2 band below the Fermi 
level and the x 2 — y 2 above. 

In the magnetic state the electron pockets become 
gapped. The peak at low energies in Bi g disappears while 
peaks at the energy of these gaps appear. Bi g samples 
the two gaps responsible for the magnetic peaks in a' xx (ui) 
and a' yy (uj). If these two gaps are different enough a two- 
peak structure should be expected in Bi g . For the values 
displayed in Fig. |9ja) these peaks arise at energies com- 
parable to that of the interband transition at T between 
the 3z 2 — r 2 band and the hole pocket, also splitted by 



magnetism, and a wide structure is observed. 

Big is less sensitive to magnetism, as also seen exper- 
imentally.^ The low energy peak from the transition be- 
tween the hole pockets at T shifts to slightly higher ener- 
gies. The one at 0.7 eV acquires a double peak structure 
and its intensity is suppressed due to the gap opening 
at the xy hole pocket at M. Spectral weight appears 
around 0.3 to 0.6 eV due to transitions between a mag- 
netic folded band a non- folded one along T — X and T — Y 
directions. 



Big samples the 3z 2 — 



— y 2 transition at T 



whose shape is strongly affected when entering in the or- 
bital differentiated regime with increasing magnetic mo- 
ment. As shown in Fig. |9^b) the spectrum changes con- 
siderably in this regime. Due to the modification of the 
x 2 —y 2 band shape, shown in Fig. [5] the 3z 2 — r 2 x 2 —y 2 
transition acquires a peak shape instead of a step one. 
Big is less affected while the spectral weight is shifted to 
higher energies. 



V. DRUDE WEIGHT ANISOTROPY 

Transport experiments in detwinned samples in the 
(tt, 0) magnetic state have shown larger resistivity in 
the ferromagneti c ^-direc tion than in the antiferromag- 
netic x-directionJ25H2Hini The anisotropy is sensitive to 
disorder and doping. It changes sign in some hole 
doped samples.^ Propo sals to explain the anisotropy 
include orbital ordering) 31 * 34 ! 37 ! integrated or at the 
Fermi surface,^ spin nematicityj^ anisotropic Fermi 
velocities,^! and anisotropic scattering from disorder in 
the anisotropic magnetic stated Fittings of the Drude 
peak performed to study whether the anisotropy orig- 
inates in the Fe rmi v elocities or in the scattering rate 
are controversial! 16 * 29 ^ The magnetic state is anisotropic. 
Even if the scattering has a strong effect on the resistiv- 
ity anisotropy, understanding the anisotropy originated 
in the reconstructed band structure is necessary. 

Within the mean-field framework above and using the 
band velocities at the Fermi surface we previously^ cal- 
culated the anisotropy of the Drude weight assuming a 
regular tetrahedron ctFe-As — 35.3°. For not very large 
intraorbital interactions the experimental anisotropy was 
found. With increasing interactions and magnetic mo- 
ments the sign of the anisotropy changed. The experi- 
mental sign of the anisotropy was concomitant with the 
smaller values of orbital ordering, discarding orbital or- 
dering as the mechanism for the observed anisotropy.^] 
We proposed that the experimental anisotropy originated 
in the topology and morphology of the Fermi surface. 
Later we argued that in the orbital differentiation regime 
the system gains kinetic energy in the ferromagnetic y- 
direction inducing a sign change in the anisotropy.^ 
Here we study the sensitivit y of the anisotropy to small 



To| shows the Drude ratio 



changes in the latticed Fig. 

D x /D y calculated using Eq. (11) for three different val- 
ues of a, the angle between the Fe-As bond and the Fe 
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FIG. 7. (Color online) Upper figures: Energy bands for the tight-binding model of Ref. [5T] in the non-magnetic (a) and (n, 0) 
antiferromagnetic state (c) for U = Jh = and U = 1.42 eV and Jh = 0.2517 respectively. Linewidths and colors reflect the 
orbital content yz=re&, 2a;=green, xy—blue, 3z 2 — r 2 =orange and x 2 — j/ 2 =black. Lower figures: Zoom of the band structures 
in the upper figures in the energy region close to the Fermi level, (b) The allowed transitions are the same as in Fig. [3] (b) . 
(d) Transitions involving a non-folded and a folded band allowed in this case but not in Fig. p](d), see text. 



plane: slightly squashed, regular and slightly elongated 
tctrahedra with aFe-As — 33.2° (left), ape-As — 35.3° 
(center) and et Fe _ As — 37.2° (right) respectively. On 
spite of the different expressions used to calculate the 
Drude weight, the anisotropy map of the regular tetra- 
hedron (middle figure) found here is remarkably similar 
to that found in our previous work. This similarity sup- 
ports the interpretation of the Drude weight anisotropy 
in terms of the Fermi surface velocity.^ 

As expected, for large interactions well into the or- 
bital differentiated regiorP^ D x /D y < 1, blue color, for 
the three values of the Fe-As angle. On the contrary, 
for small and intermediate interactions the Drude weight 
anisotropy depends strongly on the Fe-As angle. The ex- 
perimental sign of the anisotropy D x /D y > 1, in yellow 
to red colors, is largely suppressed in the squashed tetra- 
hedron case. This sign of the anisotropy is, on the other 
hand, preferred in the elongated tetrahedron case. The 
values of the Drude ratio D x /D y are also much larger 
in the elongated tetrahedron than in the regular one. 
D x /D y > 1 happens even in a relatively large area of 
the phase diagram in the orbital differentiated region be- 
fore it switches sign for larger interactions. We believe 
the values of the interactions relevant for iron pnictides 
are in the region of the phase diagram with anisotropy 
sensitive to the Fe-As angle. Such a dependence suggests 
that at least at the level of the reconstructed band struc- 



ture the resistivity anisotropy is not a robust fingerprint 
of the underlying electronic state. 



VI. DISCUSSION 

In this work we have shown the potential of using the 
velocity and Raman vertices to disentangle the orbital de- 
gree of freedom in the optical conductivity and Raman 
spectrum of multi-orbital systems. The idea follows the 
use of Raman vertices in cuprates to differentiate nodal 
and antinodal regions in k-space. The k-dependence of 
the vertices depends on the symmetry of the orbitals in- 
volved in the transition. 

We have applied this method to interacting five orbital 
models for iron superconductors denned in the one-iron 
unit cell. The optical conductivity and Raman spectra in 
the magnetic and non-magnetic states have being calcu- 
lated with the band structure of these five orbital models 
treated at the mean field level. These mean field bands do 
not include the renormalization or finite lifetimes due to 
interactions. Consequently, quantitative agreement be- 
tween the calculated spectrum and the experiments is 
not expected. On the other hand, the vertex analysis is 
valid independently of the approximation used to calcu- 
late the bands and valuable qualitative information can 
be obtained. 
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FIG. 8. (Color online) (a) a xx (uj) (top) and (b) a yv (uj) in the 
orbital differentiated region for U = 1.8 eV and Jh = 0.25(7 
(black), U — 1.6 eV, Jh = 0.25(7 and electron filling ra = 5.6 
(red) using our tight-binding model (labelled tb), and for U = 
1.53 eV and Jh = 0.25(7 and n = 6 (blue) using the Graser 
et al tight- binding modef^ (labelled tt>2). A scattering rate 
of 20 meV is used. The Drude weights are tr^(u> = 0, U = 
1.8eV,n = 6) = 1.6723, <r£*(w = 0, U = 1.6eV,n = 5.6) = 
1.1 and oif (w = 0,(7 = 1.53 eV,n = 6) = 1.05, a™(u} = 
0,(7 = 1.8eV,n = 6) = 2.5277, <t(^(w = 0, (7 = 1.6eV,n = 
5.6) = 4.9817 and a'^{u = 0,(7= 1.53 eV, n = 6) = 0.45 
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FIG. 9. (Color online) (a) Bi g (black) and B29 (red) Raman 
spectra. Low frequency Raman spectrum for U = Jh = 
in the non-magnetic state (dotted lines) and for U = 1.6 eV 
and Jh = 0.25(7 in the itinerant magnetic state (solid lines), 
(b) Main figure (inset) B\ g (-B29) Raman spectrum for U = 
Jh = (non-magnetic state), U = 1.6 eV and Jh = 0.25(7 
(itinerant magnetic state) and U = 1.8 eV and = 0.25(7 
(magnetic state, orbital differentiated regime) . In B\ g the 
shape of the spectrum changes when entering into the orbital 
differentiated region. The transition between the 3z 2 — r 2 and 
x 2 — y 2 at r is marked with an arrow. A scattering rate of 20 
meV is used. 



We have seen that interband transitions involving one 
or both hole pockets at T contribute to the optical con- 
ductivity in the far and mid-infrared. For non-negligible 
scattering rates it can be difficult to separate these tran- 
sitions from the Drude peak. These results confirm pre- 
vious calculations^ which alerted against using the ex- 
tended Drude model to an alyze the optical spectrum of 
iron superconductors!^22l a feature in the optical spec- 
trum of K-dopcd Ba-122 at frequencies 50-250 cm -1 has 
been interpreted in terms of a pseudogapP^^ precursor of 
the superconducting gap. Previous works had measured 
similar features and discussed them in terms of interband 



transitions^ or localized states.^ Our calculations show 
that the transition between the two hole pockets at T 
is optically active but remains hidden below the Drude 
peak, see also Ref.[m As it is allowed only in a very small 
region of k-space it gives a very narrow contribution to 
the optical conductivity. Experimentally it could show 
up at frequencies comparable to that of the supercon- 
ducting gap, so its presence should be considered when 
discussing the spectrum. We note that this transition 
is active in the Bi g symmetry (in the one-Fe unit cell). 
This fact could help clarify the nature of the obser ved 
50-250 cm -1 feature in the optical conductivityPSSSl 
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FIG. 10. (Color online) Drude weight anisotropy D x /D y as a function of interactions for three different Fe-As angles, 
ape-As = 33.2° as in the slightly squashed tetrahedra found in LaFeAsO (left) qfe-As = 35.3° as in a regular tetrahedra 
(center) and a^e-As = 37.2° corresponding to a slightly elongated tetrahedra (right). For small and intermediate values of the 
interaction D x /D y strongly depends on the lattice structure. 



Experiments in the magnetic state show a suppressed 
conductivity below ~ 700 cm -1 and the appearance of a 
peak around 1000 cm" 1 p^] . fl2UT4 t fl6 l fl8 | andl2Tj. Bump 
like features show up around 350 cm . In detwinned 
samples the Drude peak is anisotropic, <t' xx (uj) > cr yy (oj) 
at low energies but the anisotropy reverses at higher 
energies! 16 * 21 * 29 ^ The Raman spectra show a similar sup- 
pression and a bump in all the symmetries, and a peak 
only in B\ g symmetry in the onc-Fe unit cell (B 2g in the 
FeAs unit cell)! 3 

As discussed in Sec. |IVB| in the itinerant regime the 
magnetic peaks in a' xx (ui) and a' yy (uj) respectively sam- 
ple the gaps opened at the electron pocket at Y and X 
(folded in the hole pockets at M and T) via a transi- 
tion between two folded anticrossing bands. The dif- 
ference between peaks thus measures the electron pock- 
ets gaps anisotropy and is not a consequence of or- 
bital orderl^il The experimental sign of the anisotropy 
has bee n repro duced previously in several theoretical 
works^ M 36 * 56 * 57 ! but had not been explained in these 
terms. This anisotropy and even the possibility that any 
of these transitions is forbidden is sensitive to details of 
the underlying band structure, see Sec. |IV B| We have 
shown that the anisotropy of the Drude weight due to 
the band reconstruction in the magnetic state is also very 
sensitive to small changes in the lattice structure. 

The low energy B\ g Raman spectrum samples the elec- 
tron pockets. These become gapped in the magnetic 
state, hence the B\ g spectrum shows peaks at the en- 
ergy of the corresponding gaps. These peaks are ex- 
pected around the same energies as observed in a' xx (ui) 
and a' yy (u;). If these gaps are close enough in magnitude 
a single peak instead of two would show up in the ex- 
perimental Big spectrum. Some care is required when 
interpreting the spectrum as an interband transition be- 



tween 3z 2 



band and the hole pockets at T, active 



in Big could be close in energy. This interband transi- 
tion is also affected by magnetism, especially because it 
involves the hole pockets which stop being degenerate. 
Due to the orbital symmetry B 2g does not sample the 
magnetic gaps at the electron pockets. Therefore we do 



not expect a peak in B 2g at these frequencies. Our re- 
sults are compatible with experiments. In the magnetic 
state Big shows a peak at an energy similar to the one at 
which the peaks in <j' xx {u)) and cr' yy (uj) are observed while 
there is no peak in B 2g at this energy. 

The bump like features around 350 cm' 1 in optical 
conductivity come most probably from transitions be- 
tween a folded and a non- folded band, see also Ref. l35l 
Even if experimentally bumps appear at similar ener- 
gies in <j' xx (lo) and a' yy (uj), following the vertex analysis 
we believe that they originate in different regions of k- 
space. Low energy transitions along T — X contribute to 
a' (uj) and those along T—Y contribute to a' xx (uj) . Given 
the anisotropy along these two directions, the bumps at 
a' xx (ui) and <j' yy {uS) are not expected to show equal spec- 
tra what it is compatible with experiments.^ 

B 2g samples partially the excitations that we have pre- 
viously assigned to the bumps in <t' xx (uj) and a' yy (uj). 
Thus the presence of the bump in experiments in this 
symmetry is consistent within our expectations. On the 
other hand we do not expect these excitations to be ac- 
tive in Big. The feature observed experimentally in Bi g 
around these energies should have different origin. 

In the orbital differentiation region shown in the mag- 
netic mean field phase diagram, see Fig. [4j the spectrum 
is strongly modified with a general shift to higher en- 
ergies with no clearly identifiable feature except for the 
Big signal. For yz to become a half-filled gapped state 
at the orbital differentiation transition there is a shift 
of the yz orbital to higher energies that modifies the 
3z 2 — r 2 — > x 2 — y 2 transition at T active in Bi g . As 
a result the step feature typical of the itinerant regime 
becomes a peak. 



Finally the 3z 2 



xy interband transition at M is 



active in B 2g Raman symmetry what could help clarify 
whether any of these bands cross the Fermi level and 
complement photoemission measurements. 
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